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We consider thermodynamic and dynamic phase transitions in plaquette spin models of glasses. 
The thermodynamic transitions involve coupled (annealed) replicas of the model. We map these 
coupled-replica systems to a single replica in a magnetic field, which allows us to analyse the resulting 
phase transitions in detail. For the triangular plaquette model (TPM), we find for the coupled-replica 
system a phase transition between high- and low-overlap phases, occuring at a coupling £*(T), which 
vanishes in the low-temperature limit. Using computational path sampling techniques, we show that 
a single TPM also displays “space-time” transitions between active and inactive dynamical phases. 
These first-order dynamical transitions occur at a critical counting held Sc{T) > 0 that appears to 
vanish at zero temperature, in a manner reminiscent of the thermodynamic overlap transition. In 
order to extend the ideas to three dimensions we introduce the square pyramid model which also 
displays both overlap and activity transitions. We discuss a possible common origin of these various 
phase transitions, based on long-lived (metastable) glassy states. 


I. INTRODUCTION 


As supercooled liquids approach their glass transitions, 
one observes a very sharp increase in their viscosities and 
structural relaxation times. The physical mechanism un¬ 
derlying this slow dynamics remains controversial HHa. 
Some theories, particularly the random hrst-order transi¬ 
tion theory [5] , propose that glassy systems are approach¬ 
ing some kind of thermodynamic phase transition, with 
associated collective (slow) dynamics. The existence of 
such phase transitions can be probed by computing the 
free-energy of a pair of coupled copies (or replicas) of the 
system, and searching for a transition as a function of 
both temperature and coupling strength. These transi¬ 
tions link an equilibrium-like phase where the replicas are 
different from each other (the liquid) to one where they 
become very similar (the glass) [6]. The similarity bew- 
teen the configurations is measured by an overlap vari¬ 
able, which is the order parameter for the transition. An 
alternative approach, that of dynamical facilitation [7], 
links glassy behaviour to a dynamical “space-time” phase 
transition. This transition is explored through distribu¬ 
tions of time-integrated observables, which quantify ac¬ 
tivity in the dynamics mis]. Based on these distribu¬ 
tions, one may infer the existence of transitions between 
an active dynamical phase (the equilibrium liquid) and 
an inactive phase (the non-equilibrium glass). The order 
parameter for these transitions is the dynamical activity. 

In this work, we investigate plaquette spin models of 
glasses m, for which both overlap-fluctuations and dy¬ 
namical activity-fluctuations can be analysed, by a com¬ 
bination of analytical and computational methods. We 
concentrate on two models, whose relaxation behaviour 
is similar to that of the facilitated East model - 

their relaxation times increase faster than an Arrhenius 
law at low temperatures, but the equilibrium relaxation 
time is finite at all positive temperatures, diverging only 


as T ^0. We present evidence that these models sup¬ 
port both dynamic and thermodynamic phase transitions. 
In the thermodynamic case we consider a coupling be¬ 
tween two annealed replicas, and transitions oeeur only 
for non-zero (positive) values of the eoupling. 

We argue that these results provide a connection be¬ 
tween the (apparently quite different) ‘thermodynamic’ 
and ‘dynamic’ theories of the glass transition. This con¬ 
nection is built on the idea of metastability, which is in¬ 
trinsically connected to glassy behaviour. The formation 
of a metastable state in a finite dimensional system re¬ 
quires that small perturbations in that state do not grow: 
the system prefers to relax back into the metastable state. 
This stability to small perturbations may be described in 
terms of an interfacial cost that acts to penalise local per¬ 
turbations. Different theories ascribe different origins to 
these interfacial costs, which might be either static or 
dynamic, depending on the system of interest, and the 
kinds of fluctuation being considered. However, the ex¬ 
istence of these interfacial costs seems quite generic, and 
may be useful for rationalising different kinds of phase 
transition in these systems. 

The main results of this work are as follows. We anal¬ 
yse the triangular plaquette model (TPM) in two spa¬ 
tial dimensions [HUS], and a three-dimensional variant 
of this model, which we refer to as the square pyramid 
model (SPyM). In Section we show that two (an¬ 
nealed) coupled replicas of these systems can be mapped 
to a single replica in a magnetic field, and we derive 
some useful features of this single-replica system, which 
place constraints on the kind of phase transitions that 
can occur. Some of these results were derived in previ¬ 
ous work prHIS] . but our analysis contains several new 
insights. In Section IIP we show numerical evidence that 
the TPM in a magnetic field supports a phase transition 
in the 2d-Ising universality class. It then follows from 
the mappings in Section |n| that the coupled replicas of 
these systems also support a similar phase transition. In 
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Section we show that the TPM also supports dy¬ 
namical “space-time” phase transitions, similar to those 
in [siEiiin]. In Section 1^ we introduce the SPyM, and 
show evidence that it supports phase transitions in the 
coupled replica setting, and dynamical space-time phase 
transitions. Finally in Section we discuss the re¬ 
lationships between the thermodynamic and dynamical 
phase transitions that we have found, and we consider 
the consequences of these results for theories of the glass 
transition. 


II. PLAQUETTE MODELS, COUPLED 

REPLICAS, AND MAPPING TO SYSTEM IN 
FIELD 

A. Models 

We consider plaquette spin models defined in terms 
of classical Ising spins on regular lattices, with energy 
functions of the form 

( 1 ) 

fl 

where = ±1 with i indicating a lattice site (i = 
1,...,A), and where the interactions are in terms of 
products of spins (Ji^cFj^ • • • (Tk^ around the plaquettes 
of the lattice. See [T0l[Tll[T5] for a more general overview 
of the relevant properties of these systems. On a square 
lattice, one labels each square plaquette with an index 
/i, and {cTi^, CTj^,..., dk^} is the set of four spins on the 
vertices of plaquette ja. This construction is easily gener¬ 
alised to higher dimensions: for a cubic lattice and cubic 
“plaquettes”, each term in the energy would involve eight 
spins. This motivates us to define plaquette variables 
’ ' ’ ^k^ • 

An interesting model in this class is the TPM ID], 
where the lattice is triangular and the interactions are 
between triplets of spins in the corners of upward point¬ 
ing triangles, 

= - ^ E (TPM) (2) 

/i.=A 

Our analysis rests on a correspondence between config¬ 
urations of the spin variables di and the plaquette vari¬ 
ables r^. If we first consider rhombus-shaped systems 
whose linear size is an integer power of 2, with peri¬ 
odic boundaries, then there is a one-to-one mapping be¬ 
tween spin configurations and plaquette configurations. 
(It is clear that every spin configuration corresponds to 
a unique plaquette configuration, but the existence of a 
spin configuration corresponding to every plaquette con¬ 
figuration is less trivial [ini na na ng.) For systems of 
different sizes or with different boundary conditions, the 
correspondence is not perfectly one-to-one, but these de¬ 
viations turn out to be irrelevant in the thermodynamic 
limit. In Section |V| below, we will also discuss the SPyM, 


a three-dimensional model with the same one-to-one cor¬ 
respondence, on the body-centred cubic (bee) lattice. 

In cases where the one-to-one mapping holds exactly, 
the fully polarised state = 1 Vi is the unique ground 
state of 0- In terms of the plaquette variables, the 
ground state is = 1 V/i, and the elementary excita¬ 
tion is a “defect”, = —1. Two-body spin correla¬ 
tions vanish in these models m, although higher-order 
spin correlations are finite and allow access to a growing 
length scale at low temperatures m- Also, since there is 
a one-to-one mapping between spins and plaquettes, the 
thermodynamic properties of these models are those non¬ 
interacting binary plaquette variables [TOl [T2j [T4j [15] , or 
a free gas of ‘defective’ plaquettes (with = —1) [21]. 

However, while the thermodynamic properties of pla¬ 
quette models are trivial, their (single spin-flip) dynamics 
is not. This effect arises because flipping a single spin di 
changes the states of all of the plaquettes in which it 
participates. The plaquette dynamics is therefore “ki- 
netically constrained” [ini na na Eg, possibly leading 
to complex glassy dynamics at low temperatures. This 
is what occurs for example in the TPM whose dynam¬ 
ical properties are similar to those of the East facili¬ 
tated model [ini na ng, displaying “parabolic” super- 
Arrhenius relaxation, dynamic heterogeneity, and other 
characteristic features of the glass transition [4]. 


B. Coupled replicas 

To probe thermodynamic overlap fluctuations, we con¬ 
sider two coupled replicas of a plaquette model [ilTlIll- 
Eg • The energy function of the combined system is 

O-'') = Ej{( 7°-) + Ej{a^) - e E 

i 

where d^ and d^ are the spin configurations in the repli¬ 
cas a and b. The overlap, 

(4) 

i 

measures how similar the two copies are, and the strength 
of their coupling given by its conjugate field 5. The 
coupling Q is denoted annealed since both replicas are 
allowed to fluctuate on an equal footing. The case of 
quenehed coupling, in contrast, involves one of the repli¬ 
cas being frozen in an equilibrium configuration. Here we 
will only consider the case of annealed coupling which is 
easier to treat both analytically and numerically. Hence 
the partition function for these two coupled replicas is 

Z2(J,£) = y] (5) 

cr“,cr^ 


where the sum runs over the configurations cr^,cr^: that 
is, over all erf = ±1 and all = ±1. Here and in 
the following, we sometimes set /3 = 1 where there is no 



3 



FIG. 1. (Color online) Sketch of the duality relation in the 
TPM. The left panel shows two TPM systems, a and b, with 
coupling as in §. The right panel shows the location of the 
sites of the dual problem, again two coupled TPMs a* and 
5*. The plaquettes in the dual system bisect the coupling 
interactions in the direct system, and vice-versa, giving rise 
to the duality transformation (14). 


ambiguity [for example, the left hand side of © should 
strictly be Z 2 (/ 3 J, but we suppress the dependence on 
for simplicity]. 


C. Mapping to single system in field 


Performing the sum, we arrive at 

W2{q\J,e) = (4cosh/?J)^/2 .g-/3£;^/(9)+/3£Ei9i. (7) 

with 

PJ' = logcosh(/3J) ( 8 ) 

We recognise the exponential term in 0 as the statistical 
weight of a configuration a = g' for a single plaquette 
model with energy scale J', in a magnetic field h = e. 

To explore the consequences of this property for the 
free energy, we observe ^ 2(^5 so that 

Z 2 {J,e) = (4cosh/3J)^/^ • Zi{J',s), (9) 

where 

Zi{J,h) = E 

a 

is the partition function of a single plaquette model in a 
field h. In addition, this latter system is known to have 
an exact duality nnun] 

Zi{J,h) = {smh/3Jsmh2/3h)^/^Zi{J,h), (11) 


In [18], a mapping was derived between the free en¬ 
ergies of the coupled system 0 and a single plaquette 
model in a magnetic field. Here, we present a mapping 
between (sets of) configurations of these systems, which 
extends that analysis, as well as recovering the same map¬ 
ping between free energies. 

We introduce overlap variables qi = crfcJ^ on each site: 
our aim is to calculate the statistical weight of a partic¬ 
ular configuration of these variables. This weight is 

W 2 {q\J,e) = - <d)- (6) 


where 

= tanh(/3/i), = tanh(/3J/2). (12) 

From Eqs. the duality of the coupled plaquette 

system follows: 

Z 2 {J,e) = {smh/3Jsmh/3e)^Z2{J*,e*), (13) 

with 

= tanh(;5J/2), = tanh(/ 35 / 2 ). (14) 


We now perform the sum over the a variables. If we sum 
over (7^ first we obtain, 

W 2 {q\J, e) = ^ exp 

cr® L M 

i 

For the summation over we replace cfi^cfy ''' 
r^. Then we use the characteristic feature of the model, 
that plaquette and spin configurations are in a one-to- 
one correspondence, so we replace the sum over the erf 
with a sum over the rf. 


W 2 {q\J,s) 


^exp 


/V 

~2' 




i 


This duality is precisely the one obtained in [18] for 
the two coupled replicas of the TPM. (Note that if 
tanh^ = then tanhx = which follows from 

the definition of the tanh function, and facilitates inver¬ 
sion of these duality transforms.) 


D. Dualities and phase transitions 


The mapping from two coupled replicas to a single sys¬ 
tem in a field has useful consequences, since we may 
exploit existing results for plaquette models in mag¬ 
netic fields. The duality of this model, Eq. (11) (see 
also ini), implies a duality relation for the free energy 
Fi = - log El, 

N - - N 

Fi(h, J) ^ sinh(2/i) = Fi{h,J) F — logsinh(2/i) 

(15) 

where h and J are given in Eq. ( 12 ), and we set = 1 
as above, without loss of generality. Phase transitions 
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appear as singularities in the free energy density fi = 
limjsf^oo Fi/N. For any given J, it follows that if there 
is a single phase-transition at some h, it must happen on 
the self-dual line (h, J) = (h, J), for which 

J=—logtanhh. (16) 

On this line one has also sinh J sinh2h = 1. Such phase 
transitions were investigated in [HI HI]: the plaquette 
models considered there support a single critical point 
that occurs at some point (Jc, he) on this line, with first- 
order phase coexistence occurring on the part of the line 
with J > Jc. 

From the above mapping, these transitions correspond 
to phase transitions in the coupled replica system: the 
first-order transition line separates a state with low over¬ 
lap (small s) from one with high overlap (large s). For 
the coupled replicas, the self-dual line is 

sinh(;dJ) sinh(;d5) = 1, (17) 

This situation, where the self-dual line for the coupled- 
replica system contains a first-order transition region and 
a critical point, was proposed for the TPM in Ref. [18]. 
We present numerical evidence for this situation in Sec¬ 
tion jm] below. 


E. Other consequences of dualities and symmetries 


In this section, we explore some further consequences 
of the results derived thus far. First, we note that the 
relation 0 means that for a coupled replica system at 
parameters (J, £), the probability of a particular configu¬ 
ration of the overlap variables q is the same as the prob¬ 
ability of finding the configuration a = g for a single 
system in a field, with parameters {J'^h = s). From a 
numerical perspective, the single system in a field is much 
simpler to simulate, and the result 0 means that such a 
simulation provides direct access to all observables based 
on the overlap variables. (This result is much stronger 
than a mapping at the level of free energies.) 

Second, for a geometical interpretation of the duality 
relation (14), we refer to Figure]^ The original coupled 
system can be thought of as a lattice consisting of two 
parallel layers, a and b. The duality relation (14) may be 


interpreted as a mapping between two different two-layer 
systems, where the plaquette energy scale in one model 
determines the interlayer coupling in the other, and vice 
versa. Fig. 0b) illustrates this situation, in which the 
interlayer ‘bonds’ in the original system intersect the in¬ 
tralayer plaquettes in the dual system, and vice versa. 
This geometrical way of seeing the duality easily gener¬ 
alises to other lattices and plaquette interactions. 

Third, the duality relation for a plaquette model in a 
field can be used to analyse the behaviour of its free en¬ 
ergy in the vicinity of a (presumed) critical point. Given 
(15), it is useful to define the “singular part” of the free 


(a) Ising (b) TPM 




FIG. 2. Illustration of the relation between critical behaviour 
of the Ising model and the TPM. (a) Ising phase diagram, the 
hi = 0 axis is a symmetry line, there is a critical point indi¬ 
cated by a circle, with a first-order (phase coexistence) line for 
large J/, indicated by a solid line. Selected lines of constant 
Ji are indicated by dotted lines, (b) The corresponding situ¬ 
ation for the TPM in a field. On the solid/dashed line (16), 
the system has a discrete (Z 2 ) symmetry: a critical point and 
phase coexistence both occur on this line, as indicated. The 
dotted lines are obtained from (19) for three different values 
of /lo, and correspond to the lines of constant J/ in panel (a). 
Near the critical point, they indicate the direction of the most 
relevant renormalisation group flow. 


energy of a plaquette model in a field 

N 

Fsing(/i, J) = Fi{h, j) + ^ logsinh(2/i) (18) 


We assume that a critical point exists somewhere on the 
self-dual line, and that this critical point is in the Ising 
universality class, as is found generically [16] (see also 
below). We can then consider a renormalisation group 
flow near the critical point, which will have two relevant 
directions in the (h, J) plane. One of these relevant di¬ 
rections corresponds to the ferromagnetic coupling Jj of 
an associated Ising model; the other to a magnetic field 
hi in the Ising model. 


For the plaquette model, the relation (15) is associ¬ 


ated with a symmetry that will be spontaneously broken 
at the phase transition. Hence, the relevant direction 
that corresponds to the Ising coupling must be the one 
that preserves the self-dual symmetry: this direction lies 
along the self-dual line. However, it will be useful in the 
following to also identify the direction that corresponds 
to the Ising field hi. To accomplish this, we define a 
family of curves in the (h, J) plane that correspond to 
lines of constant J/ in the Ising system. The family of 
curves is parameterised by the points (ho, Jo) where they 
cross the self-dual line. The curves are invariant under 
the duality transformation (just as lines of constant J/ 
are invariant under the Ising symmetry transformation 
h —h.) Furthermore, the singular part of the free en¬ 
ergy, calculated along each curve, is an even function of 
h — ho (just as the Ising free energy is even in h/, for 
constant J/). 
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FIG. 3. (Color online) Simulations of the TPM in a field, (a) Distribution of the magnetisation at various values of J for state 
points on the self-dual line at system size L — 128. The bimodal distribution P(m) indicates a first-order transition, 

which disappears on reducing J. From 0, the same distributions would be obtained when considering the overlap between 
two coupled TPMs, at appropriate state points, (b) Representative configuration at phase coexistence (/3J = 2.9 and L — 128) 
showing interfaces between regions of small and large magnetisation (corresponding to regions of small and large overlap in 
the two-replica problem), (c) At our estimated critical point, (Jc = 2.634, he = 0.072) and for various system sizes, we show 
distributions of the variable x that is obtained by rescaling the order parameter M to zero mean and unit variance. The full 
line is the corresponding result for the 2d Ising model at criticality [2^, indicating that the critical point of the TPM in a field 
(and therefore of the two coupled TPMs) is in the 2d Ising universality class. 


The relevant family of curves is given by the formula 


^hoW = -logtanh(2/io - h). (19) 


The dual of any point on a line Jho{h) is easily verified 
to be (J^h) = {Jho {2ho — h) , 2/io — h) , which also satisfies 
( [l^ , confirming that the line is invariant. It then follows 
from that Fging is a symmetric function of h — ho. 
The resulting situation is shown in Fig. 

The significance of the curves ( p^ is that they give 
the direction of the most relevant RG flow at the critical 
point. The most natural order parameter for the phase 
transition is then obtained by taking a derivative of the 
free energy along these lines, for which 


d 

d^ 


d 

dh 


dJho 

dh 




We therefore define the order parameter 


= - A [i? + (V J/2)] =M -2Nd sinh J, 


( 20 ) 


( 21 ) 


where M = —{dF/dh) = is the magnetisation, 

and Nd = {d/dJ)[F + (VJ/2)] = i ^^(1 - r^) is the 
number of defective plaquettes. A key advantage of this 
order parameter is that only even cumulants of A4 may 
show singular behaviour at phase transitions: the odd 
cumulants can be shown to depend only on derivatives of 
the non-singular log sinh 2/i term in the free energy. The 
use of this order parameter also accounts automatically 
for “field-mixing” m when analysing critical properties. 


III. NUMERICAL RESULTS FOR THE TPM IN 
A FIELD 

We performed numerical simulations of the TPM in 
a field, to analyse its phase behaviour. Working always 


on the self-dual line (ITtI we use continuous time Monte 
Carlo simulations |28l I29j to sample a reweighted Boltz¬ 
mann distribution P(cr) (x where h{M) 

is a bias function. We measure the resulting distribution 
Pb(M) of the magnetisation, but we choose the function 
h{M) so that this sampled distribution does not include 
any deep minima (free energy barriers) [30]. The ‘true’ 
distribution P(M) associated with the unbiased model is 
then easily obtained as P(M) ex Pb(M)/6(M). 

Figure [^a) shows numerical evidence that for large 
J (and small /i), the self-dual line is associated with a 
first-order phase transition: see also Ref. El- The figure 
shows the distribution P(m) of the magnetisation density 
m = M/N, for three state points on the self-dual line. 
At large J the distribution is bimodal, characteristic of 
first-order coexistence. A typical configuration at these 
conditions is shown in Fig. 2(b), showing coexistence of 
low- and high-m regions, separated by sharp interfaces, 
as expected for a first-order transition. Based on smaller 
systems, a previous study speculated that phase sep¬ 
aration would not occur for the TPM in a field, but our 
results show that this does indeed occur large enough sys¬ 
tems are considered. Given the mapping 0,Fig.||b) is 
also a representative configuration of the overlap between 
two coupled replicas, for suitable (J, e). 

As J is decreased (or equivalently, temperature and 
field are increased) along the self-dual line, the bimodal¬ 
ity in P(M) becomes less pronounced and eventually 
disappears. Fig. [^a) indicates that the first-order line 
terminates at a critical point {Jc, he) with Jc ~ 2.6. 
To identify the universality class of this phase transi¬ 
tion, we performed a finite-size scaling analysis, using 
the order parameter M. defined in (21). Note that if 
all spins are up, one has M = N and = 0^ giving 
A4 = N. On the other hand, in a state with h = 0 then 
M = 0 and = 1/(1 + e*^), at low temperatures this 
gives Ai ~ —N. In general, one expects a crossover be- 














6 


TPM (2D) SPyM (3D) 



FIG. 4. (Color online) Phase diagrams of coupled plaque- 
tte models, the two-dimensional TPM (left) and the three- 
dimensional SPyM (right). The full line corresponds to a line 
of hrst-order transitions between a thermodynamic phase of 
small overlap and one of large overlap between the replicas. 
This curve is on the self-dual line (17) (dashed line). The 
hrst-order transition line ends at a critical point that is in 
the 2d Ising universality class for the TPM and the 3d Ising 
universality class for the SPyM. 


tween these two limits as h is increased from 0, with the 
crossover occuring near A4 = 0. Figure 2(c) shows the 
distribution of the order parameter M at our numerically 
estimated critical point, (Jc = 2.634, he = 0.072), for 
various system sizes L. These distributions are rescaled 
to have a mean of zero a nd a variance of unity: the vari¬ 
able X = (A4 — — {Ad)Y). The key point 

is that the shapes of the curves do not depend on the 
system size, as expected at a critical point. 

These critical distributions are universal: Fig. I^c) also 
shows the corresponding data for a two-dimensional Ising 
model at the critical point [26] . The distributions match 
well, indicating that the system is in the Ising univerality 
class. We also calculated the ratio of the susceptibilities 
X = at criticality, for the system sizes L = 128 

and L = 256. We find x{L = 2b6)/x{L = 128) = 3.36. 
Theory predicts that this ratio should scale as 
where Yh is the exponent associated with the field con¬ 
jugate to the order parameter m- For 2d-Ising univer¬ 
sality, one has Yh = 15/8 which yields a prediction of 
2^/^ ^ 3.364 for the ratio of susceptibilities. Hence, these 
results are also consistent with 2(i-Ising universality. 

Bringing together these results, we arrive at the phase 
diagram shown in Fig. We re-introduce the tempera¬ 
ture T = l3~^ as an explicit parameter and plot the phase 
diagrams as a function oiTjJ and ejJ since this is con¬ 
ventional representation in supercooled (glassy) liquids. 
The form of this phase diagram was proposed in [T8] . 
However, the results here are not consistent with the 
arguments presented in that work for the location and 
universality class of the critical point. Our finding that 
the critical point is in the Ising class is interesting, since 
this is the expected result from other general arguments 
m, and what is observed in simulation of coupled liquids 


[24l[25]. However, a key feature of the plaquette models 
is that the first-order transition line does not intersect 
the e = 0 axis except at T = 0 [T8] . 


IV. ACTIVE-INACTIVE DYNAMICAL 
TRANSITIONS IN THE TPM 


The TPM falls into a category of glassy models that are 
thermodynamically simple but where glassy behaviour 
arises because of non-trivial dynamical pathways to the 
equilibrium state at low temperature miiiiiiisiisi]. In 
fact, the dynamics of the TPM [15] is closely related to 
that of a two-dimensional East model [IIHlS]- Many ki- 
netically constrained models, including the East model, 
display dynamical phase transitions - phase transitions 
in the space of trajectories - between a phase with a 
high dynamical activity, A, and one with low dynamical 
activity [8]. Eor these lattice models, the activity K is 
defined as the total number of configuration changes in 
a trajectory [32] [33]. By coupling a field s to the dy¬ 
namical activity one can define the so-called 5-ensemble 
(also known as the exponentially biased or tilted ensem¬ 
ble) where the probability of obtaining a trajectory X^- 
of total time extension r, is re-weighted by its activity 
00132] 


Ps{Xr) 


e-^^Po{Xr) 

Zsir) 


( 22 ) 


Here Ps{X^) is the probability of the trajectory X^ in 
the 5-ensemble, Po{Xr) is the unbiased probability of this 
trajectory (the one generated by the actual dynamics of 
the system), and Zs{r) is the moment generating func¬ 
tion of the activity A, in effect a partition sum for tra¬ 
jectories. Eor long times Zs(r) takes on a large-deviation 
form UMm. 


Z,(r) (23) 


where 0{s) plays the role of a dynamical free-energy, and 
is the scaled cumulant generating function of the activity. 
The mean time-intensive activity of trajectories in the 5- 
ensemble. 


k{s) = {K)Jt, 


(24) 


where the average is w.r.t. the ensemble (22), can thus 


be obtained from k{s) = Dynamical - or “space- 

time” - phase transitions manifest as singularities in 0{s) 
0111321135]. We also define the susceptibilty 


X(s) = = T ^{ 5 K^) 


(25) 


Like the East model, the dynamical relaxation of the 
TPM is hierarchical, due to energy barriers to relaxation 
that are logarithmic in the linear size of relaxing regions 
[12]. This in turn leads to a “parabolic” [36] super- 
Arrhenius law for the typical relaxation time in equi¬ 
librium as a function of temperature. Given the similar 
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dynamical properties of the TPM and the East model 
dynamics, a natural question is whether the TPM also 
displays active-inactive “space-time” transitions. 

In order to answer this question numerically we make 
use of transition path sampling (TPS) [37] to efficiently 
sample trajectories in the 5-ensemble. For the purposes 
of numerical efficiency we exploit the fact that different 
trajectory ensembles can be defined by fixing different dy¬ 
namical quantities, such that for long trajectories these 
distinct ensembles become equivalent jSS] [39|. This is 
analogous to what occurs with standard configuration en¬ 
sembles in equilibrium statistical mechanics in the large 
size limit. We consider in particular the x-ensemble intro¬ 
duced in [39], that is, the ensemble of trajectories with 
fixed number of configuration changes, i.e. activity, K 
but where the overall trajectory time extension r fluc¬ 
tuates. In contrast, the 5-ensemble above is one where 
trajectories are of fixed overall time but their activity 
fluctuates. For large K and r these ensembles can be 
shown to be equivalent [39l|40]. For the case of the TPM 
the x-ensemble is particularly efficient to simulate (see 
[39] and m for details) and the functions k{s) and x(5) 
can be recovered from the ensemble equivalence. 

Figure shows the results of the 5-ensemble analysis 
of the TPM. It shows the average activity for a system of 
size = 8 X 8 as a function of 5, at temperature T = 0.5 
(note this is the TPM in the absence of field). As the 
length of the trajectories is increased the change in {K)s 
becomes more pronounced, as seen in the corresponding 
susceptibilities. This is indicative of a first-order transi¬ 
tion at some Sc ^ 0. Similar size scaling is observed by 
changing system size, as shown in the insets. The depen¬ 
dence of this transition on the temperature is discussed 
in Sec. jV C| below, together with similar results for the 
(three-dimensional) SPyM. 


V. OVERLAP AND ACTIVITY TRANSITIONS 
IN A THREE-DIMENSIONAL PLAQUETTE 
MODEL 

In order to explore whether the static and dynamical 
transitions found above for the TPM are present in di¬ 
mensions other than two it is of interest to generalise 
the TPM to higher dimensions. One of the reasons is 
that if one wishes to consider plaquette models to study 
“quenched” coupled replicas [a Ham] or “random pin¬ 
ning” [42H44| the distinction between two and three di¬ 
mensions may be very significant, due to the inability 
of two-dimensional systems in random fields to support 
first-order static transitions [45]. Here we introduce a 
three-dimensional model that is similar to the TPM. 


A. Model 

The model we consider is defined on a three- 
dimensional body-centred cubic (bcc) lattice. The “pla- 



A 

FIG. 5. (Color online) Average activity k{s), and the asso¬ 
ciated susceptibility x(s) in the TPM at T = 0.5. The main 
panels show data for system size L = 8. These results were 
obtained by the a:-ensemble method (see text), using trajecto¬ 
ries with fixed numbers of events A, as shown. On increasing 
the trajectory length, the crossover from active to inactive 
behaviour becomes increasingly sharp and the susceptibility 
peak increases. The behaviour for smaller systems (L = 4) 
is shown in the insets, with both quantities normalised by 
the system size. In the absence of a phase transition, one 
expects both k{s)/L^ and x(s)/L‘^ to be independent of L, 
so the sharper crossover at L = 8 is again consistent with an 
underlying phase transition. 

quettes” are upward-pointing square-based pyramids, 
each containing five spins. Considering the standard bcc 
unit cell, one such pyramid is formed by the spin at the 
centre of the cube together with the four spins at the cor¬ 
ners of the lower face: see Fig.|^ We call this model the 
square pyramid plaquette model, or SPyM. Its energy 
function reads, 

Ej{a) = ~^ai^aj^ak^ai^am^ (SPyM) (26) 

/X 

where /i runs over all the pyramidal plaquettes on the 
lattice, and the location of the five interacting spins 
shown in Fig. [bl This is “model 1” of [16]. 

Just like the TPM, the SPyM has a one-to-one cor¬ 
respondence between spin and plaquette configurations. 
An alternative model [inmni may be defined on a face- 
centred cubic lattice, in which the plaquettes are tetra¬ 
hedral pyramids (“model 2” of m)- However, in this 
case each interaction involves four spins so the system 
has a global spin-flip symmetry, and the spin-plaquette 
correspondence is not exact in finite (periodic) systems. 
However, these deviations from the one-to-one correspon- 



















FIG. 6. (a) The SPyM consists of spins (grey circles) on 

the sites of a BCC lattice, which interact in quintuplets at 
the vertices of upward pointing square pyramids. One such 
pyramid is indicated; the central spin also participates in four 
other upward pointing pyramids whose apexes are the four 
spins on the upper face of the cube. 


dence are irrelevant in the thermodynamic limit. 

Returning to the SPyM, we explicitly demonstrate the 
one-to-one correspondence between spins and plaquettes, 
by a general method that applies also to the TPM. We 
choose as basis vectors for the lattice di = (1,0,0), 
a 2 = (0,1,0) and as = (—1, —1, v/2)/2. We focus 
on systems whose sites are at Idi md 2 — nds with 
l^m^n G {0,1, 2,..., L — 1}, with periodic boundaries 
(so for example sites with n = L — 1 are neighbours of 
those with n = 0). We indicate the location of the ja-th 
pyramid by the position of the spin at the apex. The 
plaquette variable for fi = (i, j, k) is then 

^ ^(i —l,jf —1,/c—1) • 

Following the same reasoning as in m we can invert 
this relation in terms of a “Pascal pyramid”: the idea is 
to demonstrate that introducing a single defect into the 
system corresponds to flipping a particular set of spins. 

Starting from the ground state, we demonstrate the 
procedure by introducing a single defect at the origin: 
this affects those spins in upper layers which lie on the 
sites of an inverted Pascal pyramid (or fractal pyramid). 
Assuming that the central spin in Fig. is at the ori¬ 
gin, we flip that spin, which introduces a defect in the 
pyramid below it. In order to avoid any other defects, 
we also flip the four spins on the top face of the cube 
shown in Fig. which ensures that there are no defects 
in any of the pyramids pointing upward from the origin. 
Iterating this procedure for all other layers, the final spin 
configuration is 




^ ^ mod 2 


(28) 


where (^) = are combinatorial numbers, and 

0 < i, j < k (all other spins have a = 1). Given periodic 
boundary conditions, this procedure determines all spins 
in the system: on setting the final layer of spins, it may 
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T In t 

FIG. 7. (Golor online) Relaxation of the energy of the SPyM 
at low temperature starting from a random configuration (the 
system size is L = 16). The curve shows the character¬ 
istic plateaus indicative of hierarchical relaxation, as in the 
East model and the TPM. Inset: average relaxation time as 
a function of inverse temperature, showing super-Arrhenius 
behaviour. 


be that defects in the final layer are unavoidable. How¬ 
ever, for systems whose linear size L is a power of 2, it 
is easily shown that this procedure produces a final state 
with exactly one defect. Now observe that for any spin 
configuration, flipping the set of spins for which a = —1 
in (28) inverts the state of the plaquette variable just be¬ 
low the origin, leaving all other plaquette variables con¬ 
stant. Similarly, to flip the state of any other plaquette, 
one applies a spatial translation to set of spins for which 
a = —1 in (28), and flips all the spins within this trans¬ 
lated set. Hence, by repeatedly applying this procedure, 
one can generate a spin configuration that corresponds to 
any given configuration of the plaquette variables. This 
is sufficient for a one-to-one correspondence between spin 
and plaquette configurations, since we know already that 
the total number of spin and plaquette configurations are 
the same, and that every spin configuration corresponds 
to exactly one plaquette configuration. 


This one-to-one correspondence between spin and de¬ 
fect configurations means that the thermodynamics of 
the SPyM is that of a free binary gas of plaquettes. Fur¬ 
thermore, the relaxational dynamics is similar to that of 
a (three-dimensional) East model. Figure shows the 
decay of the energy at low temperatures starting from 
a T = oo configuration. We see the characteristic hi¬ 
erarchical decay of both the East model and the TPM: 
the energy decays in steps with characteristic time scales 
Tn = with n = 0,1,2,.... These steps become ap¬ 
parent when plotting as a function of the rescaled time 
variable (T/J) log t [46]. The inset to Fig. shows that 
the equilibrium relaxation time of the SPyM is super- 
Arrhenius, as in the East model and the TPM. 
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(b) TPM (c) SPyM 



In Fig. ^b,c), we show how the crossover in activity 
varies with the temperature T/ J, for both the TPM and 
the SPyM. Simple estimates [13 SH] indicate that if the 
inactive state is metastable and relaxes to equilibrium 
via some kind of nucleation process with rate 7nuc per 
unit volume, then s* Jnuc/^k, where 6k is the activity 
difference (per unit space-time) between the active and 
inactive states m- We attribute the existence of the 
transition in this model to a stable inactive state with 
almost no defects. We expect that the rate for relaxation 
back to equilibrium is a strongly decreasing function of 
temperature, which is consistent with the increasing s* 
as T/J increases. (The activity difference Ak between 
active and inactive states increases with T/J but this 
dependence is much weaker than that of the relaxation 
rate.) 


VI. DISCUSSION 

A. Connection of phase transitions to long-lived 
metastable states 


FIG. 8. (Color online) (a) Activity as a function of s in the 
SPyM for system size L = 4, various trajectory lengths, and 
T/J = 0.65. The inset shows the corresponding susceptibility, 
(b-c) s-ensemble phase diagrams for the TPM and SPyM. The 
full curves are an estimate of the transition point from the 
simulations. The dashed lines are extrapolations in the low 
temperature regime inaccessible to numerics. 


B. Phase transition in (annealed) coupled replicas 

The SPyM possesses the exact duality described in 
Sect. II.C. In particular, the properties of the SPyM in a 
field were studied in Ref. m (“model 1” of that paper), 
where it was found that on the self-dual line there is a 
first-order transition between phases of small and large 
magnetisation terminating at a critical point in the 3D 
Ising class. From those results we can directly infer the 
phase diagram of the two coupled SPyMs via the map¬ 
ping of Sect. II.C. The result is shown in Fig. This 
phase diagram is similar to that of the TPM, except that 
the range of phase coexistence is larger and the critical 
point occurs at higher temperature. 


C. Evidence for a dynamical (space-time) phase 
transition 

As well as the phase transition for coupled replicas 
in the SPyM, we also present evidence for a space-time 
phase transition, similar to that shown for the TPM in 
Fig-i The results for the SPyM are shown in Fig.[^ for 
temperature T/J = 0.65 and linear system size T = 4. 
There is good evidence for a sharp transition at s = s* > 
0, as found in the TPM. 


The phase behaviour shown in Figs. [^ and reveals 
striking similarities between thermodynamic transitions 
(for coupled replicas) and dynamic transitions (based on 
dynamical activity). We now argue that these transitions 
are connected to the existence of long-lived metastable 
states, which are intrinsically linked to the glassy be¬ 
haviour in these systems. 

Consider first transitions for annealed replicas. If we 
work at the phase coexistence point, but within the high- 
overlap phase, the system occupies low-energy states. 
For 5 = 0, these states would be metastable. This 
metastablility means that when localised low-overlap re¬ 
gions are generated by thermal fluctuations, these re¬ 
gions tend to shrink, just as small fluctuations tend to 
shrink within classical nucleation theory. Escape from 
the metastable state requires a collective process that 
operates on some finite length scale. As 5 is increased 
from zero, this length scale increases, as does the associ¬ 
ated free energy barrier: both diverge at the coexistence 
point where the high-overlap phase becomes stable. 

The situation for dynamical phase transitions is sim¬ 
ilar, except that one should think of trajectories of the 
system as (J+l)-dimensional objects that exist in space- 
time. If one works at the dynamical phase coexistence 
point (some 5 = 5* > 0) then inactive trajectories dom¬ 
inate the s-ensemble. During these trajectories, the sys¬ 
tem remains localised in low-energy metastable states, 
with small thermal fluctuations of the activity, associ¬ 
ated with space-time “bubbles” [7]- If the trajectory 
length r is less than the time required for escape from 
the metastable state, similar inactive trajectories can be 
generated with unbiased (s = 0) dynamics, by taking ini¬ 
tial conditions from low-energy metastable states. Trans¬ 
formation of such a trajectory into a typical equilibrium 
trajectory involves the introduction of an active space- 
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time bubble, which subsequently grows to macroscopic 
size. The connection with metastability arises because if 
one introduces a small active bubble within an inactive 
trajectory, one expects to incur a cost in probability (if 
this were not the case, the state would not be metastable 
since it would readily relax back to equilibrium). As in 
the case of overlap fluctuations, the critical bubble size 
and the probability barrier increase as s is increased from 
zero, diverging at the coexistence point. 

We argue that this analogy between phase coexistence 
phenomena induced by s- and e-fields provides a qualita¬ 
tive explanation of the similarity between Figs. and 
in that both are linked to the existence of metastable 
states that can be observed in unbiased (s = 0 = e) sys¬ 
tems. The relevant metastable states have low energy: 
both the inactive state of Fig. and the high-overlap 
state of Fig. have much lower defect concentrations 
than the equilbrium average value c ~ . For large 

s, the system minimises its propensity for dynamical ac¬ 
tivity by removing defects, so that Nd/N ^ 0; for large 
coupling 5 it is easy to show that N^/N ^ 
the system has an effective temperature 1/(2^) [49]. 

Evidence for phase coexistence induced by s and 
5-fields have both been presented in atomistic mod¬ 
els laEn. By contrast, in kinetically constrained mod¬ 
els, dynamical phase coexistence can be induced by the 
5-field |8| but there is no such transition as a function 
of 5. Metastable states do exist in these systems [48] . 
so phase coexistence at positive s may be expected, as 
argued above. However, this metastability does not lead 
to a phase coexistence for any positive 5 - the reason is 
that the “dynamical” metastability that occurs in KCMs 
appears because thermal noise only generates a certain 
subset of the overlap fluctuations that are possible within 
the metastable state. If all fluctuations in the overlap 
are possible, one finds that localised low-overlap regions 
grow easily in these systems. This means that there is no 
barrier between active and inactive states in the coupled- 
replica construction. It is only within the dynamical 5- 
ensemble that the metastability becomes apparent, due 
the restriction on the kinds of fluctuation that are gen¬ 
erated by the thermal noise. 

Finally, it is important to note that the inactive and 
high-overlap states in Figs. and are structurally dis¬ 
tinct from their equilbrium counterparts. This is quite 
different from ensembles with a quenched coupling be¬ 
tween replicas m, where the structures of high- and 
low-overlap states should be statistically (almost) indis¬ 
tinguishable. 


B. Connection between multiple coupled replicas 
and biased-activity ensembles 



FIG. 9. Three dimensional stack of coupled two-dimensional 
TPMs. 


posed of many replicas of the d-dimensional system ar¬ 
ranged parallel to each other along the extra dimension, 
see Fig. Such system of n coupled replicas has an 
energy 

E„{a\a\ ...,) = Ej{a^) + Ej{<j^) + ■■■ 

E +■•■)• (29) 

i 

Using the methods of Refs. nanHiisi] it is easy to prove 
that the partition sum of the (d+l)-dimensional problem 
also has an exact duality: 

Zn{J,e) = (sinh^Jsinh/3e)^"Z„(J*,e*), (30) 


where we have assumed periodic boundary conditions in 
the transverse direction, and (J, 5 ) and (J*,5*) are re¬ 
lated again by (14). Similar results were found in [52| 


for other classes of plaquette models. Given this duality 
we expect the phenomenology of this many-replica sys¬ 
tem to be similar to that of two replicas, except that any 
phase transitions should be in the (d + 1)-dimensional 
Ising universality class (assuming that both longitudinal 
and transverse dimensions are taken to infinity in the 
thermodynamic limit). 

The partition sum Zn has a natural transfer matrix 
representation in the transverse direction, = Tr(T’^), 
with 


T = cosh^(/3J/2)e^^" 


1 + tanh(/3J/2)cr|)^cr|^ • • • 


(l + e-'^Vf), 


(31) 


where are Pauli matrices. This in turn can be 

related to the generator of (imaginary time) quantum 
evolution in the usual manner [52| when the trans¬ 
verse coupling is large and the longitudinal one small 
jg-2/3e, tanh(/3J/2) ^ 1], so that Z^ oc exp(—tid), with 




i 




cr\ 


ku ’ 


(32) 


where 


The connection between Figs. [^ and can be further 
motivated through a generalisation of the coupled two- 
replica system discussed above. Given a plaquette model 
in dimension d, consider the associated d+1 system com- 


6t h = e St g = /dJ/2, t = nSt. (33) 

The Hamiltonian ( j^ generates dynamics in the trans¬ 
verse direction. While it is not derived from a stochas¬ 
tic operator it has the basic features of the generator 
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la sa for the dynamical ensemble defined in an 

off-diagonal part (the terms) that perform configura¬ 
tion changes and a diagonal part (the plaquette terms 
associated to the escape rate. The parameter s in the 5- 
ensemble operator controls the relative strength of the 
diagonal and off-diagonal terms HEg, in analogy with 
the balance between h and g in ( [3^ . Furthermore, the 
duality (30) implies a duality h ^ in (32), with the pos¬ 
sibility of a dynamical transition at that self-dual point 
g = h. This connection between a static transition in the 
1-dimensional problem ( [2^ (itself closely connected to 
the static transition in the two-replica plaquette system) 
and a dynamical transition in the d-dimensional system 
(32) provide another rationalisation of the similarities be¬ 
tween Figs. [^and[^ 


C. Outlook 

Plaquette spin models have several features that make 
them attractive for studies of the glass transition. As 
we have shown here, exact results can be derived, which 
guide numerical studies of phase behaviour and many- 
body correlations. The models are also computationally 
much less demanding than atomistic models of super¬ 
cooled liquids, so that (for example) finite size scaling 
over a large range of system sizes can be performed, to 
analyse phase transitions. The equilibrium relaxation of 
the models follows a dynamical facilitation scenario, in 
which point defects play a central role. However, there 


are strong many-body correlations, and the statistics of 
overlap fluctuations are rich and complex, as anticipated 
in the theory of Franz and Parisi [6]. In this sense, the 
models provide a bridge between different theories. In¬ 
deed, as argued in m, one might describe plaquette 
models by a modified form of RFOT, but with two im¬ 
portant caveats: (i) the analogue of the Kauzmann tran¬ 
sition occurs at zero temperature in these models (ii) the 
interfacial cost associated with growing droplets of a new 
state within a typical equilibrium state scales logarithmi¬ 
cally in the droplet size (not as a power law as anticipated 
by RFOT). 

Looking forward, we hope that further work on pla¬ 
quette models (particularly in d = 3) will show to what 
extent mean-field [6] and RFOT ideas can be modified to 
apply in this setting. We can imagine that the apparently 
different physical pictures envisaged by thermodynamic 
and dynamical theories of the glass transition [4] might 
both be applicable in these models. In that case, it is 
not clear whether some new results would be required to 
discriminate between the theories, or whether they might 
in fact offer complementary descriptions of the same phe¬ 
nomena. 
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